Towards a Quantum Fluid Mechanical Theory of Turbulence 
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Recent studies of turbulence in superfluid Helium indicate that turbulence in quantum fluids obeys 
a Kolmogorov scaling law. Such a law was previously attributed to classical solutions of the Navier- 
Stokes equations of motion. It is suggested that turbulence in all fluids is due to quantum fluid 
mechanical effects. Employing a field theoretical view of the fluid flow velocity, vorticity appears as 
quantum filamentary strings. This in turn leads directly to the Kolmogorov critical indices for the 
case of fully developed turbulence. 
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I. INTRODUCTION 

It is often thought that fluid mechanical turbulence 
is the last of the great unsolved problems of classical 
physics*]]]. Even with extensive numerical solutions of 
the classical Navier-Stokes fluid mechanical equations of 
motion, it has not been possible to prove simple exper- 
imental rules, which have been formulated for turbulent 
fluid flows @. Our purpose is to suggest that turbulence 
is inherently a quantum mechanical problem which can- 
not be solved by classical methods. 

To see what is involved in a quantum mechanical pic- 
ture, let us consider the role of quantum mechanics in de- 
scribing vorticity and turbulence. The classical Reynolds 
number describing a turbulent states is defined as 



K 



VC _ VCp 
v rj 



(1) 



wherein V and C represent, respectively, the velocity and 
length scales of the flow, v and r\ represent, respectively, 
kinematic and dynamic viscosity of the fluid and p is 
the fluid mass density. Turbulent flows are normally de- 
scribed by high Reynolds number. The number of quan- 
tum vortex filaments characterizing a flow as described 
by Feynman[l[ reads 



n F = , 

K 



(2) 



wherein the quantized vortex circulation k for a fluid con- 
taining molecules each of mass M and total atomic mass 
number A is given by 



2-kU _ 3.96 x 1(T 3 cm 2 
M ~ A sec 



(3) 



For the formation of a single quantum vortex in an oth- 
erwise circulation free fluid flow, the Feynman number 
Rp ~ 1. The ratio 

7c K 

denotes the number of quantized filaments within a tur- 
bulent tube normally described employing classical vor- 
ticity. For water H20[i| with A = 18 and at room tem- 
perature and pressure, 



K 



v 'water 
^"water 



45. 



(5) 



Vortex tubes in turbulent water flows may then be safely 
considered to be constructed from quantized vortices. 

In Sec|II]the quantized theory of fluid mechanics will 
be reviewed § and formulated as a Lagrangian quan- 
tum field theory. Quantized vortex filaments will then 
be discussed in Sec lIIIl In SeclIVI the general theory of 
quantum fluid quantization will be written in the Cleb- 
sch representation p. In Sec [V] the statistical correlation 
functions of velocity and vorticity are defined and dis- 
cussed. The classical Kolmogorov model@,Si for fully 
developed turbulence will then be considered wherein the 
critical index for velocity fluctuations will be derived. In 
Sec lVIl the Kolmogorov scaling exponents will be de- 
rived under the assumption that quantized vortex fila- 
ments have the geometry characteristic of a self avoiding 
random walk. When a path is embedded in D spatial 
dimensions, the fractal dimension of a random walk path 
is two, while the fractal dimension of a self-avoiding ran- 
dom walk path is 

, D + 2 
d = — - — , 



5/3 for D = 3. 



(6) 
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The fractal dimension d = 5/3 of a random quantized 
vortex line uniquely determines the Kolmogorov scaling 
exponents for fully developed turbulent quantum fluid 
dynamics. In the concluding Sec lVIIl the very similar 
nature of turbulence in the quantum superfluid phase of 
4 He2 and 3 He—B and in ordinary fluids usually regarded 
as classical will be discussed. It will be argued that tur- 
bulence in fluid flows normally regarded as classical are 
in reality inherently quantum mechanical in nature. 



II. QUANTUM FLUID MECHANICS 

Let M a be the mass of the a th molecule at position r a 
in a fluid. The mass density p and the mass current J 
operators may be defined 

p(r) = ^l«(r-r a ), 

a 

J(r) = i£(p B *(r-r B )+*(r-r a )p B ), (7) 

a 

wherein the momentum p Q conjugate to the position r a 
obeys 



[p a ,r b ] = -ihlS ab 



(8) 



The mass density and the mass current thereby obey the 
equal time commutation relations 



[J(r) )/9 (r')]=p(r)grad5(r-r') 



(9) 



One may formally define the fluid velocity operator v 
according to 

J(r) = ~(p(r)v(r)+v(r)p(r)) (10) 

in which case the velocity field does not commute with 
the mass density field; i.e. 



[v(r),p(r')] = grad 5(r — r') 



(11) 



Note the consequences of the Jacobi identities which fol- 
low from Eq. (|ll[) in the form 



[ Vl (r)>,(r'),P(r")]] = 0, 
[p(r"),K(r'),^(r)]] = 0. 



(12) 



The mass density commutes with the velocity component 
commutators. By taking the curl of both sides of Eq. (TlTj) . 
it follows that the vorticity components also commute 
with the mas density, 



fl(r) = curlv(r), 
[n(r),p(rO]=0. 



The commutation relations between different compo- 
nents of the velocity read 

l -p{v) [« i (r),« J -(r')] = ^[«i(r) J « J -(r')]p(0 

l -p{v) [v i (r),v j (r')}=n ij (r)S(r-r'), 

Qij(r) = e ijk tt k (r). (14) 

Let us now consider the quantum vorticity and fluid cir- 
culation in more detail. 



III. VORTEX FILAMENTS 

For a given velocity potential function 3>(r), consider 
the unitary operator 



U = exp ( ^ / p(r)$(r)d 3 r 



Acting on a wave function \l/ of N fluid molecules, 
'iM N 



C/* = exp — £$(r a ) 



(15) 



(16) 



The unitary operator U comutes with the mass density, 

rfp(r)U = p(r), (17) 

but not with the velocity, 

[/ f v(r)J7 = v(r) + grad$(r). (18) 

Ea.ljlS]) follows directly from Eas.([TT]) and fTS]). Let us 
now consider the circulation of the velocity around a 
closed curve C; i.e. 



r(C) = d> v • dr 



(19) 



Eqs. fTS)) and (fT^)) imply the circulation operator tranfor- 
mation law 



U^T(C)U = r(C) + <j> grad$ • dr. 



(20) 



The last integral on the right hand side of Eq. (j20|) need 
not vanish if the closed curve C resides in a multiply 
connected region of nonzero mean fluid density. 

For example, suppose a vortex filament with a finite 
core which is empty of fluid molecules. Let C completely 
surround the core. From the single valued nature of quan- 
tum mechanical wave functions, in particular the velocity 
potential phase in 



(13) 



N 

n 



,iM*(r )/ft 



(21) 
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one must have the quantized circulation 

2nh 



grad$ • dr 



c 



M 



N(C) , 



N(C) = 0,±1,±2,±3,±4. 



(22) 



The central problem of quantum fluid mechanical tur- 
bulence concerns the nature of the motions of filaments 
with a quantized circulation in units of k = (2nh~/M). 
Let us now turn to a formal Lagrangian representation 
of a field theoretical formulation of fluid mechanics. 



IV. CLEBSCH QUANTIZATION 

Let us at first begin with the classical fluid mechanical 
action principle 10] for adiabatic flows. Later the action 
will be employed to quantize the theory, The results are 
consistent with the Landau quantization approach dis- 
cussed in the above Sec|n] From a classical fluid me- 
chanical viewpoint, the velocity field may be written in 
the form 



v = grad$ + Agrad/i 

wherein the three component velocity fields (w 



(23) 



.) 

are replaced by three scalar fields <E>, A and p. The vor- 
ticity may then be computed as 



fi = curlv = gradA x grad/i. 



(24) 



Vortex lines may be pictured as the intersection between 
two surfaces where one of the surfaces has spatially con- 
stant A and the other surface has spatially constant p. In 
adiabatic flows, one expects that the vortex lines move 
with the fluid. This expectation may be implemented by 
allowing three classical conservation laws; i.e. 



dp 
dt 



djpX) 

dt 
d(pp) 

dt 



div(pv) = 0, 
div((pA)v) = 0, 
div((p/i)v) = 0. 



(25) 



In terms of the total time derivative 
d d 

- = -+v.grad, 

we have the conservation laws in the more simple form 



(26) 



dp 
~dt 



— pdivv 
dX 

dp 
~dt 



0, 
0. 



(27) 



Thus, each of the fluid particles and their A and p surface 
coordinate assignments are conserved. The spatially con- 
stant surfaces of both A and p move with the fluid and 



so do the vortex lines which are the intersection between 
these surfaces. 

Eqs. (f2"Tl) are an expression of Kelvin's circulation 
theorem[ll| for adiabatic flows. In particular, the cir- 
culation around a closed boundary curve C = <9E of a 
surface E may be written 



no 



■ dr = j> 



Xdp) 



(28) 



wherein Eq. (f2"3")) has been invoked. If E lies within a 
simply connected fluid region of non-zero classical fluid 
density, then Stokes theorem implies 



r(C) = i Xdp = 

JC=dT, 



dX A dp. 



(29) 



If C, E, A, and p all move with the adiabatic flow, then 
the circulation T(C) is uniform in time in accordance 
with Kelvin's theorem. 

Let u(p, s) representing the fluid energy per unit mass 
as a function of the mass density p and the entropy per 
unit mass s. The classical fluid action may be written as 



.4 



Ldt 



Cd 3 i 



dt. 



with an adiabatic lagrangian density 

C = p(^-u( P ,s)]+C'. 



(30) 



(31) 



In the above Eq. (|3"Tj) . the first two terms on the right 
hand side represent the usual Galilean invariant "kinetic 
energy minus potential energy" form while the conserva- 
tion law constraint lagrangian density may be written 



dp 
~dt 



jodivv 



(32) 



The velocity potential $ appears as a Lagrange multi- 
plier assuring local conservation of mass, A appears as a 
Lagrange multiplier assuring local conservation of p and 
local conservation of A is the result of the classical vari- 
ational action equations of motion; i.e. SA = 0. Let us 
consider this in more detail. 

Functionally differentiating the action Eq. (|30|) with re- 
spect to the velocity yields Eq. (|2"3"]) in the form, 



SA 



= p (v — grad$ — Agrad/i) = 0. 



(33) 



The conservation law Eqs. ([27|) follow from the variational 
action principle equations 



SA 



dp 



SA 

Jx 



pdivv = 0, 
dX 



dp 
} ~di 



0. 



(34) 
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Still viewing adiabatic fluid mechanics as a classical field 
theory, one may compute the conjugate momenta to the 
fields of interest employing the general field theoretical 
relation 11^ = (SA/Sip) where ip = (dcp/dt). In particu- 
lar, two conjugate field-momenta may be computed via 



n = SA = dC =<p 
p Sp dp 

TT SA 9C X 



(35) 



The energy of a classical fluid flow in the lagrangian for- 
malism is 



E = 



DC DC . 

p a-p + ^-^ dr ' (36) 



which after parts integration yields 

E = J £d 3 r wherein £ = ^p\v\ 2 + pu(p, s). (37) 

The quantized fluid mechanical model may now be for- 
mulated as follows: 

(i) There are four fields of interest p, <E>, A and p,. Either 
field of the pair (p, $) commutes with either field of the 
pair (/i,A). From Eqs. (f3"5|) , we have the two conjugate 
field commutation relations pairs 



(38) 



-[$(r),p(r')]=5(r-r'), 

i r i \ \( iw S ( r ~ r ') 
-[),(,), A(r)]=-^-. 

The fluid velocity field operator is computed via Eq. 
employing the operator ordering 

v(r) = grad$(r) + ^{K r )> gradju(r)} (39) 

where {A, B} = AB + BA and with implicit point split- 
ting for field multiplications at the same spatial point. 
From Eqs. ((38|) and (|39|) . one may recover the Landau 
equal time commutation Eqs. tfTTj) . P^)l and (fl"4"|) . 

(ii) The Hamiltonian follows from Eq. ([3T]) to be 

I 



H 



Hd 3 r wherein H = ^v- (lp) ■ v + pu(p, s). (40) 



Note the operator ordering between the non-commuting 
field operators p and v. 

(iii) The fluid temperature T and pressure P may be 
computed from the thermodynamic law 



du = Tds H — -dp. 



P 

7- 



In terms of the enthalpy per unit mass, 



d(pu) 
dp 



P 

= u H , 

P 



(41) 



(42) 



we have the thermodynamic law 

dw = Tds + -dP. 

P 

For adiabatic flows ds = 0; i.e. 

p gradiu = gradP (adiabatic). 



(43) 
(44) 



(iv) The equation of motion for the mass density may be 
derived employing 



p(r) = -[ff,p(r)]=-divJ(r) 



(45) 



wherein the current density operator J is defined in 
Eq . (fTU]) . The equation of motion for the velocity field, 

v(r) = i[JT,v(r)], (46) 

has the quantum operator form 

v(r) + grad(^v(r) • v(r) + w(p(r),s)) 

+ -(n(r) x v(r) - v(r) x n(r)) =0, (47) 

with the vorticity l~i = curlv. Eq. (|4T|) is simply the quan- 
tum field theoretical version of Newton's law, 



p^- = — gradP (classical), 



(48) 



in thinly disguised form. In Eq. (|48p we have invoked 
Eqs. (|2"6")) and (|44p. Finally, the equation of motion for 
vorticity follows by taking the curl of Eq. (|47|) ; i.e. 

ft(r) = l [IT, «(■■)], 

O(r) + -curl(n(r) x v(r) - v(r) x n(r)) = 0. (49) 

This completes the local lagrangian formulation of 
quantum fluid mechanics for adiabatic flows. The quan- 
tum field theory has been well defined modulo the usual 
list of suspect mathematical rigor problems endemic to 
strongly non-linear quantum field theories describing an 
infinite number of degrees of freedom. 



V. KOLMOGOROV CORRELATIONS 

Consider the velocity fields of fully developed turbu- 
lence in a frame of reference in which the mean drift ve- 
locity are zero. The statistical correlations of the velocity 
field for isotropic incompressible turbulence (diw = 0) 
may be described by 

G(r-r') = i({v(r),v(r')}) 
d 3 k 



G(r-r')= J S(k> 



2 

ik-(r— r 



(2tt) 3 



kk 



S(k) = (l-p-]S(k), 



(50) 
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wherein {a, b} = (ab + ba) and 

G(|r-r'|) =trG(r-r') 
G(|r-r'|) = i(v(r).v(r')+v(r')-v(r)) 

sin(fer) 



G(r) = - / k 2 S(k) 
i" Jo 



kr 



dk. (51) 



The Komlmogorov energy spectral distribution £{k)dk 
is defined as the portion of the kinetic energy per unit 
mass in the wavenumber width dk; i.e. Eqs. (f50")) and (|5Tj) 
imply 

~k 2 S{k)~ 



£{k) = 

O 

G{r) = 2 / £{k) 



2tt 2 
sin(fcr) 



kr 



dk. 



(52) 



The mean squared velocity of the isotropic turbulence is 
thereby 



£(k)dk. 



(53) 



Vorticity correlations may be discussed in a similar man- 



ner 



F(r-r') = i({fi(r),fi(r')}) 



F(r 



Sn(k)e 



ik-(r— r ) 



kk 



d 3 k 



(2tt) 3 



S n (k) = \1 - S n (k), 
Sn(k) = k 2 S(k). 



(54) 



If we allow 



Fflr-r'l) = trF(r-r'), 
F(|r-r'|) = i(n(r)-n(r')+n(r')-n(r)>, 

F(r) = \T k 2 S n (k) \ s J^ll] dk: (55) 
it Jo I Kr 

then the spectral theorems hold true; 

poo 

rF(r)=2 k£(k) sin(kr)dk, 
Jo 

1 f°° 

k£(k) = — I rF)r)sm{kr)dr. (56) 
I" Jo 

The mean squared turbulent vorticity may also be com- 
puted as 



/Winn = 2 



k 2 £(k)dk. 



(57) 



The classical viscous heating rate per unit mass from 
the turbulent mass eddy currents may be taken to be 



1P> 2 

e = = vui , 



(58) 



wherein r\ is the fluid viscosity. Kolmogorov^ assumed 
that the energy spectral function depends only on the 
wave number k and the dissipation per unit mass e, i.e. 



£ = F{k,e), 



(59) 



and then considered the physical dimensions of the quan- 
tities involved. These dimensions (in cgs units) are 



dim[£] = 



cm 

sec 2 



dim[e] = , 
sec 3 

dim[fc] = — . 



(60) 



The only functional form in Eq. (|59p with the correct 
physical dimensions of Eq. (|60p is given by 



£(k)=C 



,2/3 



£5/3 



(61) 



where C is a dimensionless constant presumably of order 
unity [l2|]. The correlation between turbulent velocities 
at two different points may be defined as 

. 9 (r-r') = (|v(r)-v(r')| 2 ), 
g(r) = 2(G(0) - G(r)), 
sin(fcr 



9{r) 



£{k) 



kr 



dk, 



(62) 



wherein Eq. (|52[) has been invoked. Employing Eqs. (|6T 
and (|6"2"j) . one finds 



C" = 4C 



1 



sin/3 



dp, 



g(r) = C\erf/\ 
wherein the gamma function, 



T(z) 



, ds 
s 



(63) 



(64) 



has been invoked. Under the assumption that g = Q (r, e) 
alone, the proportionality in Eq. ([63|) . i.e. gir) oc (er) 2 / 3 
follows from the analysis of physical dimensions. From 
Eqs. (p)6")) and (|6"Tj) follows the Kolmogorov correlation 
scaling function for vorticity 



e = r G) c = 5 c ' 

. e 2 / 3 



(65) 



Let us now see how the Kolmogorov scaling Eq. ([63)1 may 
arise from quantum vortex filaments. 
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VI. QUANTUM VORTEX FILAMENTS 



which now reads 



The scaling laws in Eqs. fBlj) and (|63|) appear to 
be valid for fully developed turbulence in superfluid 
Helium [3, 0, EH, Ell as well as in liquids normally re- 
garded as classical in nature. Let us see why this may be 
so from the viewpoint of quantum fluid mechanical vortex 
filaments. From Eq. (|22|) . it follows that the circulation 
around a single quantum fluid vortex obeys 



dr 



2irh 



(66) 



The magnitude of fluid vorticity is related to the number 
of quantum vortex per unit area A^ 1 via fl — (k/A). 
If we examine just one quantum vortex filament and let 
T be the unit tangent vector to the quantized filament, 
then 



n 



A' 



(67) 



Following this single quantum vortex filament, let s de- 
note the arc length alongh the filament path so that 



T(s) 



dr(s) 
ds 



(68) 



Along a single quantum vortex filament, let us consider 
the correlation function 

H( Sl ~s 2 ) = \ <T( Sl ) • T(s 2 ) + T(s 2 ) • T( Sl )) ; (69) 

From Eqs. ([68|) and |69|) it follows that the mean square 
distance between two points on the fillament separated 
by a filament length L is given by 



R 1 



/ / H( Sl 



S2)dsids 2 - (70) 



If the filament path is of fractal dimension d, then the 
mean square radius varies with L according to 



L 
I 



2/d 



(71) 



wherein £ is a length scale along the filament core re- 
quired to achieve an appreciable bending. Evidently such 
a length should be large on the scale of the diameter of 
the filament core; i.e. £ ^> VA. In order that Eqs. ([70|) 
and (I71| hold true, one may write 



H(L) 



2(c2-l)/d 



(72) 



Employing Eqs.jnHD, ([1)7]), & and ©> one nnds for 
fractal dimension d 



f,W-(S)V)-!(£ 



r 


t 




L 



2(d-l)/d 



(73) 



2(d-l) 



(74) 



in virtue of Eqs.{7T]) and ([73j) . Eq.(74]), describing vortic- 
ity correlations by employing quantized vortex filaments 
of fractal dimension d, is the central result of this section. 

If the fractal dimension of the quanized vortex is that 
of a self avoiding random walk, as in Eq.([6]), then F(r) 
is given by 



F(r) 



-f-V 

2 \AJ 



4/3 



(75) 



which exhibits the Kolmogorov scaling as in Eq. ([55 



VII. CONCLUSION 



A quantized theory of fluid mechanics was reviewed 
and formulated as a Lagrangian quantum field theory. 
Quantized vortex filaments were discussed. A theory of 
the quantized fluid mechanics was written in the Cleb- 
sch representation. Finally, a model of quantized vortices 
was given wherein the Kolmogorov scaling law was de- 
rived. 

Recent studies of turbulence in superfluid Helium in- 
dicate that turbulence in quantum fluids obeys a Kol- 
mogorov scaling law. In this paper, it is suggested that 
turbulence in all fluids is due to quantum fluid mechanical 
effects. There have been found to be many differences (Tsl. 
Il9l . [20j between turbulence in superfluid He and other flu- 
ids, yet more striking, are the similarities plj. It has been 
found experimentally that superfluid He^ and more re- 
cently liquid helium 3 He — B exhibit the Kolmogorov 
scaling law found in other fluids, an unexpected result. 
The fact that both "classical" and quantum fluids mani- 
fest the same scaling laws, suggests that all fluids need to 
be treated using a quantum mechanical formalism when 
turbulence takes place. 

Turbulent flows occur when the Reynold's number be- 
comes large, i.e., TZ » 1. In this regime, the kinematic 
viscosity is small. Although turbulent flows are consid- 
ered highly dissipative[22|], they last a long time, even af- 
ter an external energy source ceases to exist. A vanishing 
kinematic viscosity implies that the ratio TZf /"R> — vjn 
is relatively small. The implication is that there are few 
quantum vortices in a classical bundle. The quantum vor- 
tices themselves obey the Kolmogorov 5/3 scaling law. 
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